Nanoscale magnetometry using a single spin system in diamond 



O 



R. S. Said 1 , D. W. Berry 2 , J. Twamley 1 
1 Centre for Engineered Quantum Systems, Department of Physics & Astronomy, 
Faculty of Science, Macquarie University, Sydney, NSW 2109, Australia. 
2 Institute for Quantum Computing and Department of Physics & Astronomy, 
University of Waterloo, Ontario N2L 3G1, Canada. 

We propose a protocol to estimate magnetic fields using a single nitrogen- vacancy (N-V) center 
in diamond, where the estimate precision scales inversely with time, SB ~ 1/T, rather than the 
square-root of time, SB ~ 1/ y/T. The method is based on converting the task of magnetometry into 
phase estimation, performing quantum phase estimation on a single N-V nuclear spin using either 
adaptive or nonadaptive feedback control, and the recently demonstrated capability to perform 
single-shot readout within the N-V [P. Neumann et al, Science 329, 542 (2010)]. We present 
numerical simulations to show that our method provides an estimate whose precision scales close 
to ~ 1/T (T ~ the total estimation time), and moreover will give an unambiguous estimate of 
the static magnetic field experienced by the N-V. By combining this protocol with recent proposals 
for scanning magnetometry using an N-V, our protocol will provide a significant decrease in signal 
acquisition time while providing an unambiguous spatial map of the magnetic field. 
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I. INTRODUCTION 

A highly sensitive magnetic field sensor, that can op- 
erate at room temperature and has atomic spatial reso- 
lution may revolutionize many nanotcchnologics, for in- 
stance in medical and biological technologies, advanced 
material sciences, spintronics, and quantum comput- 
ing. Toward the practical realization of such a device, 
nanoscale magnetometry experiments in solids have been 
realized using single nitrogen-vacancy (N-V) centers in 
diamond.—^ These experiments detect very weak mag- 
netic fields, ~ 3 nT at kilohertz frequencies, and can 
locate a nearby electronic spin with a spatial resolution 
of ~ 5 nm by utilizing electron spin dynamics of a defect 
center in a diamond nano-crystal, as illustrated in Fig. 
[TJ This defect center possesses remarkable properties for 
magnetic sensing: It can be individually addressed and 
optically polarized and measured, and it maintains spin 
coherence at room temperature for considerable periods 
of time i^i^ However, the current magnetometry precision 
is limited by standard statistical fluctuations, namely the 
shot-noise limit £ In this limit the precision of the mag- 
netic field estimate scales as SB ~ 1/ y/T, where T is the 
total time needed to acquire the estimate. 

This scaling is because T/t independent measurements 
are made over a short time r. This yields an uncer- 
tainty in the magnetic field o£ SB ss (h/gfiB){l/VrT)- 
In principle, if one were to use a measurement over the 
entire time interval T, then one would have a measure- 
ment with uncertainty scaling as SB ~ 1/T. This is 
the best precision possible for a measurement over this 
time interval, and is equivalent to the Hciscnberg limit for 
phase measurement^ There are two problems preventing 
measurements with precision scaling as 1/T. The first 
is spin-spin relaxation; performing a single measurement 
beyond the dephasing time T2 does not yield an increase 
in precision. The second is that performing a measure- 
ment over a longer time may result in ambiguities. That 
is, the magnetic field causes the spin to rotate more than 



once, and the number of rotations cannot be determined 
from the measurement. Experimental advances have in- 
creased T 2 , and there are proposals to extend T 2 to the 
order of a second^ 

In this work, we address the second problem, and 
present a method to achieve Hciscnbcrg-likc scaling of the 
precision for measurement times smaller than T 2 while 
eliminating any ambiguities in the estimation. This will 
allow faster acquisition of a magnetic field map for a 
preset precision. In summary, we adapt a more gener- 
alized quantum phase estimation algorithm (gQPEA)^ 
to instead estimate the phase generated by an unknown 
Z— rotation of the Bloch sphere of a qubit (atomic two 
level system), via Ramsey intcrfcromctry. The gQPEA 
phase estimation algorithm was initially developed using 
the framework of optical intcrfcromctry to estimate an 
unknown phase acquired when a photon passes through 
a static phase shifter, and it has been experimentally 
demonstrated using linear optical methodsi&S The proto- 
col can make use of either adaptive^^ or nonadaptive^ 
controls to yield unambiguous estimates of the phase with 
a precision which scales inversely with the overall mea- 
surement time, i.e., Hciscnberg- like scaling. Our protocol 
makes use of single-shot measurements (SSM) of the spin 
of the atomic system, and we generalize to the case when 
the visibility of such measurements may be significantly 
below 100%. 

We numerically simulate our magnetometry protocols 
under both ideal and realistic measurement conditions, 
taking into account atomic decay and dephasing. In the 
ideal case, when we have perfect SSM visibility, we pre- 
dict that the precision of the magnetic field estimate has 
better scaling than the shot-noise limit; i.e., the precision 
scales as SB ~ 1/T 73 , 1/2 < /? < 1. We describe this 
type of scaling as sub-shot-noise scaling. Furthermore, 
we analyze the performance of these protocols when the 
SSM visibility can be quite low, and surprisingly find 
that sub-shot-noisc scaling in the estimate precision is 
still possible. 



We begin by reviewing Ramsey interferomctry in a 
two-level system to estimate an unknown phase rotation. 
Section III adapts the gQPEA to atomic systems to op- 
erate within the Ramsey interferometry cycle. Section 
IV presents the results of numerical simulations for the 
adaptive and nonadaptive protocols with varying single- 
shot measurement visibilities. 

We surprisingly discover that, for reduced SSM visi- 
bility, the nonadaptive method performs better than the 
adaptive method. Furthermore, we find that the non- 
adaptive method can achieve a scaling close to 5B ~ 1/T 
even when the SMM visibility is substantially reduced. 
Hence, we focus our attention on optimizing the non- 
adaptive protocol to work as efficiently as possible in ex- 
perimentally realistic conditions where one may not have 
perfect single-shot measurement capability. We finally 
discuss our conclusions and possible future directions. 



II. SPIN INTERFEROMETRY 
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Before reviewing Ramsey spin interferometry, we clar- 
ify some important definitions used throughout our pa- 
per. The resource required to obtain an estimate of the 
magnetic field to a preset precision is essentially the to- 
tal interaction time between the atomic probe and the 
unknown magnetic field. We denote this total interac- 
tion time by T. For shot-noise scaling the uncertainty of 
the estimate scales as ~ 1/ ^/T, sub-shot- noise scaling is 
~ l/T@ for 1/2 < (3 < I, and Hcisenberg-like scaling is 
1/T. 

We now consider a two-level spin system which expe- 
riences a static unknown magnetic field B z . Through a 
pulse sequence, which is analogous to Ramsey atomic in- 
terferomctry, we can convert the problem of estimating 
this unknown magnetic field strength into the problem 
of estimating the unknown phase cf> oc X g B z t acquired by 
the spin during a time t, where X g is the gyromagnetic 
ratio of the spin, and t is predetermined accruing time 
between the spin and the magnetic field. Our protocol 
will estimate the value of the accrued phase 4>. 

We will choose the two interfcrometric basis states 
{|0), |1)}, whose superpositions will precess in the mag- 
netic field, to be the two nuclear hyperfine states \m s = 
0, rnj = ±1/2), as shown by the green transition in Fig. 
QJe). We choose these for two reasons: First, as recently 
demonstrated^ one can perform single-shot measure- 
ments of the nuclear spin state associated with these 
basis states, and second, these nuclear spin states pos- 
sess very long coherence times. The hyperfine split- 
ting between the two states \m s = 0,mj = ±1/2), 
arises from the coupling between the electron spin of 
the 15 N-V center (S = 1) in its ground state, and the 
15 N nuclear spin (I = 1/2). However, another choice 
might be to choose the shorter lived electronic spin states 
{\m s = 0), \m s = ±1)} of the N-V center^ as depicted 
by the blue transitions in Fig. [He). 

Our protocol requires one to perform many single-shot 



FIG. 1: (Color online) (a) A single spin of an N-V center 
(blue point), attached to the tip of an atomic force micro- 
scope probe, detects a weak magnetic field from another spin 
(red point). The 'blue' probe's dynamics are controlled by 
microwave wires (two black stripes), (b) Structure of the cen- 
ter, consisting of a substitutional nitrogen impurity N (blue 
point) adjacent to a carbon vacancy V (white point), tetra- 
hedrally coordinated to each other and three nearby carbon 
atoms C (red points).— The impurity can be a 14 N(7 = 1) or 
1 N(7 = 1/2) atom, (c) Electron spin levels of the center. The 
ground state 3 A is a triplet (S = 1) with a 2.88 GHz zero-field 
splitting between m s = and m s — ±1. Standard Ramsey 
interferometry performed in the electron spin is the basis for 
the operation of many magnetometry protocols.—^ (d) Hy- 
perfine level splitting due to a coupling to the 15 N nuclear 
spirit in the absence of nuclear Zeeman effect, (e) Degen- 
eracy (m s = 0, mi — ±1/2) can be lifted to allow optically 
detected nuclear spin resonance. Red, blue, and green arrows 
represent optical, microwave, and radio- frequency transitions, 
respectively. The qubit proposed to probe the weak magnetic 
field is the 'green' transition. 

Ramsey measurements using a range of time durations. 
We summarize each sequence as follows: The spin is ini- 
tialized to a state |0) = |m s = 0,m/ = —1/2), which 
can be pictured as a vector along the +Z direction in 
the Bloch sphere. This initialization can be carried out 
in the N-V system by optical cycling, near resonance on 
the optical transition 3 A — 3 E. We then apply a 7r/2- 
pulse on resonance with the |0) — |1) transition, and this 
produces a pure superposition state, written in density 
matrix form as 

po = ^(|0)(0|+i|0)(l|-*|l)(0| + |l>(l|). (1) 

That is, the initial vector is rotated to the equator of 
the Bloch sphere, and that rotation is assumed to be im- 
plemented in a very short time. This rotation can be 
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realized by applying a radio-frequency field if we use the 
nuclear spin to encode the interferometric basis states, or 
by using a microwave field in the case of electronic spin 
states. The spin system then undergoes free precession 
due to the effect of the unknown external static magnetic 
held B z . For some predetermined period of time t, the 
spin system accrues a relative dynamical phase (f>, pro- 
portional to the magnetic field strength, and rotates the 
Bloch vector around the Z axis in the X — Y plane. 

During this precession we have the Hamiltonian Hf = 
X g B z a z , and including decay and dephasing, we can de- 
scribe the evolution via a master equation in the Lindblad 
formic, as 



Pit) 



(2) 



where C(A,B) = ABA* - \ iftk,B^, a z = 

|1)(1|, and <j_ = |0)(1|. We denote the decay time as T\. 
Setting p(t = 0) = po, and h — 1 for simplicity, we obtain 
a solution to ([2j as 
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where A = 2i\ g B z — 1/T%. We then apply another ir/2- 
pulse (tt/2 rotation around the X axis), mapping the 
phase information into population information. We then 
perform a single-shot measurement of the population to 
get a single detection result. For clarity we call such a 
single-shot measurement a 'click', in contrast to the over- 
all measurement of magnetic field, which includes many 
of these individual 'clicks'. The probability of measuring 
the system in the state |0) or |1) (+1 or -1 click along 
the Z axis), conditioned on the unknown field B z , is 



p(±i|£y 



i 



l±e 



-t/T 2 



cos i2X g B z t) 



(4) 



For a small predetermined accruing time t = r, where 
t « ?2, and repeating the standard prepare-evolve- 
measure procedure M times, the total amount of time 
resource expended is T = Mr. This naive procedure 
provides a standard deviation of the estimate that scales 
as£ SB Z ~ 1/ \/T. This is the shot- noise limit for mea- 
surement precision. In what follows, we will show that 
by adapting the gQPEA protocol,^ which involves ac- 
quisition of 'clicks' over varying time durations, we can 
obtain an unambiguous magnetic field estimate with pre- 
cision that scales better than the shot-noise limit. 



III. ESTIMATION METHOD 

To provide clear insight into the estimation method, 
we recast (@| into 



P(±l|^) 
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where the system phase <f> = 2X g B z t, is the unknown and 
constant real quantity we would like to estimate. The 
obtained estimate itself, after the estimation procedures, 
is denoted by <f> to differentiate it from the actual un- 
derlying system phase. Following Ref£, we expand the 
probability distribution for the system phase after the 
m-th single-shot measurement as 



(0) 



after 



and can compute the expectation value for cxp(iq 
the m-th click as (e^) m = J P m (^)e^d0. 

If there is no initial knowledge about the phase, the 
initial probability distribution P m= o(^) is flat. After 
the first detection, m = 1, we have one bit of informa- 
tion about the system phase that can be used to cal- 
culate the three nonzero coefficients ^={o,±i},m=i- Af- 
ter m detections arc made, a vector of detection re- 
sults u m = {ui, U2-, ■ ■ ■ , u m } is obtained, where each 
Ui(i = 1, ...,m) = ±1. Using Bayes' rule, the condi- 
tional probability distribution for the system phase given 
the next detection result is 



P MOm+l) <X P (±l\<f>) P (<l>\ 



(7) 



with P{4>\u m ) = P m (<fi). It gives an update formula for 
the unnormalized coefficients of the probability distribu- 
tions P m (4>) H> P m+ i(4>) as 



, (8) 



where the normalized coefficient is = 

frj,m+l/27T&o,m- 

The estimate of the system phase <f> after a total of m = 
M detections is then obtained by taking the argument of 
the expectation value of e 4 ^, 



arg(e^) M 



arg 



P i4>\u m=M ) 



arg (b- 



l,M) 



(9) 



Averaging over a number of trials (or samples) S, we 
obtain an estimate of the Holevo variance Vh—, corre- 
sponding to the square of the uncertainty in the phase 
estimate, as 



(|27r&_! 



1. 



(10) 



The case when SB ~ 1/s/T, (usual shot-noise limit) can 
be recast as VjjT ~ C (a constant), while the limit where 
SB ~ l/T can be recast as V H T - 1/T. Below we will 
identify sub-shot-noise scaling where V H T ~ 1 > 

/?>0. 

Rather than using standard Ramsey interferometry 
for phase estimation we make use of another protocol, 
originally devised in an optical setting, which uses mul- 
tiple accumulations of the unknown phase to achieve 
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near Hciscnbcrg scaling of the uncertainty in the phase 
estimate The multiple accumulations correspond to 
a prolonged accruing time, which provides more informa- 
tion about the phase. 

Applying this multi-accumulation concept to our prob- 
lem, we modify the expression of (4) by replacing t = r 
with t = Tfc = 2 k r: 



P{±i\4>) 



l±e" Tfc/T2 cosf2 i 
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where now <p = 2X g B z r. Prolonging the accruing time in- 
creases the rotations of the state vector in the A— Y plane 
and consequently introduces a modulo 2n ambiguity into 
the estimate. We eliminate this ambiguity through ad- 
justment of the parameter k, and the application of an 
additional effective phase shift $ to the system. This 
can be achieved via the active application of an external 
magnetic field or passively, through simply performing 
the second Ramsey 7r/2 rotation about a rotated axis 
on the A — Y plane of the Bloch sphere instead of the 
X— axis. This additional phase shift can be chosen adap- 
tively (as feedback) or nonadaptively (as a predetermined 
phase increment), such that 

P(±l\<f>, *) = i [1 ± V(r k ) cos (2 fc - $)] , (12) 

where V(rk) — exp (— 2 A Y/T 2 ) is the estimation visibil- 
ity. Given a conditional detection probability we can 
determine the Fisher information?^ 



J~6 



1 



d(f> 



4 fe sin 2 (2 fe - $) 
exp (2 1+fe r/T 2 ) - cos 2 {2 k cj) - $) : 



(13) 



where £ is a click measurement result. 

The Fisher information essentially represents the 
amount of information about cj) contained in the mea- 
surement results £. Maximal information about the sys- 
tem phase can be extracted when we choose values of k 
that maximize (|13[) . for a fixed ratio of t/T 2 . Further 
investigation of (fT3|) allows one to show that no useful 
information regarding the system phase can be inferred 
when T = 2 k T 3> T 2 . This is numerically confirmed by 
our simulations as presented below. The overall measure- 
ment time can be greater than T2, but the uncertainty 
scales as 1/vT for T > T- 2 . Beyond T 2 we can only hope 
to achieve an uncertainty which scales, at best, like the 
shot-noise limit. 

In the following, we first describe a protocol for sub- 
shot-noise magnetometry which uses adaptive feedback, 
where the value of the control phase $ in (fT2|) is de- 
termined based on previous detection results. Following 
this we describe a method which is perhaps more suitable 
for experiments, where the control phase $ is regularly 
incremented without any dependence on the previous de- 
tection results. 



A. Adaptive scheme 

For the adaptive control scheme, the control phase $ in 
(fT2| is updated after each detection based on the previ- 
ous detection results. Ideally, this control phase is altered 
in such a way so as to minimize the variance of the fi- 
nal phase estimate. However, there is no known method 
to determine such a phase analytically, and numerical 
methods are computationally intensive. Instead, we use 
an adaptive update scheme^ to minimize the variance of 
the phase estimate after the next detection, based on the 
information obtained from the previous detections. 

In order to minimize the variance after the next detec- 
tion, one must maximize the quantity .M— , 



1 

2^ 



E 

«m =0,1 



(14) 



The adaptive control phase that maximizes M. is one o£i£ 

$0 = arg (ba* — c*a) , 

(15) 



$± = arg 



ci 



where 



ci = (a*c) 2 -(a&*) 2 + 4(|6| 2 -|c| 2 )&*c, 

c 2 = -2ilm(a 2 b*c*), (16) 

and a, b, and c are functions of the &j, m _i expansion co- 
efficients for the probability distribution, a = b-\ t m —i, 
b = i>_i_2*,m-i^( r fc), c = b -i+2*, m -iV(-Tk)- Finally, we 
arrive at the update formula for bj <m n- bj <m +i, based 
on the acquisition of the (to + l)-th click measurement 
Um+i) and using (|12p and the control phase <I>o,± that 
maximizes Ai, to obtain 



m+1 



+ U m+ iV(Tk) 



>j-2 k ,m t 



-i*n 



L) j+2 k ,m 
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(17) 



For clarity, we enumerate the individual steps in the 
estimation procedure using the adaptive control phase as 
follows: 

1. All parameters are set: the fundamental ac- 
cruing time r, the maximum number of multi- 
accumulations 2 K , the maximum number of detec- 
tions M, the number of trials S, and the dephasing 
time T 2 . 

2. For the s-th trial, the protocol starts using the ini- 
tial k — K and the initial additional phase $. 

3. The spin state is initialized to |0). 

4. The 7r/2-rotation around the A-axis is applied by 
the external pulse (on-resonance with the |0) — |1) 
transition) to create the coherent superposition 
state. 
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5. The Bloch vector undergoes free evolution in the 
X — Y plane for t% = 2 k r, due to the non-varying 
unknown magnetic field and the additional phase 
<F 

6. Another 7r/2-rotation around the X-axis is applied 
to transform the phase information into the spin 
population. Instead of the active application of the 
phase shift $ in step 5, one can simply rotate the 
axis of this second Ramsey rotation. 

7. A single-shot measurement is performed to provide 
the result £ = ±1. 

8. The phase estimate 4> is inferred and the additional 
phase <J> is updated based on the measurement re- 
sult u m . 

9. Steps 3-8 are repeated M times. 

10. Steps 2-9 are repeated for k = K.K — 1,K — 
2,...,0. 

11. Steps 2-10 are repeated S times. In the end, one 
has S number of final estimates required to obtain 
a numerical estimate of the Holevo variance Vh ■ 

To complete a single trial (steps 1-9), the protocol re- 
quires a total time resource T a f = M x (2 K+1 — 1)t. 
Below we will use T 2 /r = 10 3 , and taking^ T 2 ~ 2 ms, 
we have a fundamental accruing time of r = 2 ^ts. For 
simplicity, we will from now on make reference to dimcn- 
sionless time quantities, which have been rescaled by the 
fundamental accruing time as T = T/t, such that 

f af = M (2 K+1 - 1) . (18) 

B. Non-adaptive scheme 

Updating the control phase $ without any dependence 
on the previous detection results requires that number of 
detections vary as a function of k£ From Refi^, we choose 
a simple linear function for the number of detections, 

M{K,k) = M K + F(K-k), (19) 

where Mk is an initial number of detections at k = K, 
and F is a positive integer. The number of detections 
increases as k is decreased. Hence, the total time resource 
required to complete the protocol with the nonadaptive 
scheme is 

f na = M (2 K+1 -1)+F (2 K+1 -2-K). (20) 

We determine the control phase $ at the m-th click mea- 
surement as 

$ m = $ m _l+7T/2, (21) 

and use the same updating formula for the b's as (|T7|). 
Furthermore, the individual steps of this scheme are sim- 
ilar to those of the adaptive one, except for the following 
steps: 
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FIG. 2: (Color online) Logarithmic graphs of the total time 
resource T a f versus the quantity VnTaj of the phase esti- 
mate. We plot the shot-noise scaling as the black dashed line 
and ~ 1/T (Heisenberg-like) scaling as the black solid line. 
The numerical result of our protocol with adaptive feedback 
is plotted with a red dotted-solid line. For T a f < T2, we 
achieve Heisenberg-like scaling of the estimate precision. We 
set the number of samples S = 1000, the number of single- 
shot detections M — 6, K = 1 — > 20, and the dephasing time 
f 2 = T 2 /r = 10 3 . For f af > 10 3 , the method fails and we 
can, at best, achieve shot-noise scaling through increasing the 
number of detections M at a fixed K. 

1. All parameters are set: the fundamental ac- 
cruing time r, the maximum number of multi- 
accumulations 2 K . the initial number of detections 
Mk, the number of trials S, and the dephasing time 
T 2 - 

2. For the s-th trial, the protocol starts using the ini- 
tial values of k = K , Mk , and $. 

8. The phase estimate <p is inferred and the control 
phase $ is updated following (|21l) . 

9. Steps 3-8 arc repeated for M (K, k) = M K +F(K- 
k) times. 

In contrast to the adaptive scheme, where one has to 
run the protocol from the maximum k = K to the lowest 
k = 0, the nonadaptive scheme can be operated with- 
out adhering to the order. This flexibility could possibly 
simplify an experimental realization. 

IV. SIMULATION 

We numerically simulated the above protocol for the 
adaptive scheme, in the case of perfect SSM visibility, by 
setting the number of trials S = 1000 and the number of 
single-shot measurements M = 6. The simulations were 
performed for values of K from 1 to 20. Meanwhile, the 
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FIG. 3: (Color online) Logarithmic graphs of T a f(T na ) ver- 
sus VnT a f (VnTna) when the detection visibility is imperfect 
fd < 1, in the case of (a) the adaptive scheme, and (b) the 
nonadaptive scheme. The red (blue) dotted-solid lines set 
f d = 95% (f d = 85%), and set K = 1 -> 20, S = 5000, and 
T- 2 /t = 10 3 for both simulations. We choose F = 2 here for 
the nonadaptive method. 

ratio of dephasing time to t is chosen to be T2/T = 10 3 . 
A random phase guess is used as a prior estimate to ini- 
tiate the protocol. The simulation result is presented 
in Fig. [2] where we plot the simulated quantity VjjT a f 
against T a f on logarithmic scales. As can be seen from 
Fig. [2j the precision of the simulated phase estimate has 
Heisenberg-like scaling for T a j < 10 3 . If one proceeds 
with the estimation protocol for T a f > 10 3 , the precision 
will be worse than the shot-noise scaling (not shown in 
the figure). Therefore, to maintain the shot- noise scal- 
ing for T a f > 10 3 , one needs to stop the protocol at a 
certain value of K that gives T ss T2, and subsequently 
increase M only. It is important to note that for these 
measurements, no prior knowledge of the actual phase is 
needed. 

To consider a more realistic condition where a perfect 
single-shot measurement is unattainable, we quantify the 
imperfection both by a detection contrast and a detection 
visibility. The detection contrast of SSM is quantified by 
two parameters: a maximum and a minimum of the de- 
tection probability to find the spin state in an eigenstate 
of the detection basis. Here, we denote such maximum 
and minimum probability to get the +1 click by f a and 



fi. Using these parameters, we straightforwardly gener- 
alize (dU) into 

p(±l\ < f ) ^j a j i ) = k±Ji±l^Ay, (22) 

where y — V(rk) cos(2 fe — <f>), giving the updating for- 
mula for the probability distribution coefficients as 

bj, m +l = (f a + fi) bj tm + r " +1 (f a - fi) V(r k ) 

x ('', 2 . .,„< + 6 i+2fc , m e itE> '™) . (23) 

For the more symmetric case where f a + fi = 1, we set 
f a — fi = fd, and simplify (|2"2")l such that 

P(±l\<f>, $, = X - (1 ± f d V(Tk) cos (2*0 - <&)) , (24) 
giving the updating formula 

bj,m+l = bj^ m + U m+ ifdV(Tk) 

x ( bj-2K m e~ l ' 5 ' m + b^k^e 1 *™ ^ ^ 

The detection visibility is included in our simulations 
to compare the robustness of both adaptive and non- 
adaptive schemes against the SSM imperfection. We first 
choose parameters to optimize the nonadaptive protocol 
for a given detection visibility (i.e. the symmetric case 
where f a + ft = 1). In the next subsection, we optimize 
the nonadaptive protocol to improve its robustness for 
imperfect detection contrast (the nonsymmetric case). 

Using (|24[) and ([231) . simulations were performed for 
a range of nonunit detection visibilities. As shown in 
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FIG. 5: (Color online) Logarithmic graphs of T„ a against 
VnTna, obtained from the optimized nonadaptive protocol 
for different detection contrasts and for some combinations 
of {M K , F} = {8, 8} (a), {12, 8} (b), {16, 8} (c), {8, 12} (d), 
{12, 12} (e) and {16, 12} (f). The detection contrast param- 
eters used are /, = 5% and f a — 55% (blue), 65% (green), 
75% (red), and 85% (cyan). The black solid line represents 
Heisenberg-like scaling, while the shot-noise scaling is shown 
by the black dashed line. The simulations use 10000 samples, 
T 2 /r = 10 3 , K = 1 -»• 8. 



Fig. [3j it was found that the nonadaptive scheme is more 
robust against imperfect detection visibility. The results 
for the nonadaptive method with a range of visibilities are 
presented in Fig. [4) indicating that a detection visibility 
of at least 90% is required by the nonadaptive scheme 
to obtain Heisenberg-like scaling. While performing the 
simulations to obtain both Figs. [3]and|ll we did not limit 
the value of K, in order to show that the variance has 
worse scaling than the shot-noise limit for T a f^ na > T2. 



Optimized protocol with nonadaptive scheme 

The nonadaptive protocol can be optimized to work 
with nonideal detection contrast. We optimize the pro- 
tocol by numerical simulations based on (|22|) and (f23|) to 
find realizable combinations of the number of detections 
M, and the multiplier F to maintain the sub-shot-noise 
scaling. This is useful to provide some optimized param- 



eters for a possible experiment where the measurement 
efficiency is often not symmetric, as quantified by the 
detection contrast {f a ,fi}- 

We present the simulation results in Fig. [5j where the 
protocol is optimized for some values of Mk = {8, 12, 16} 
and F = {8, 12}, at a range of visibility parameters f a = 
{0.55,0.65,0.75,0.85}, where U = 0.05 is fixed. 

From the figure, the overall results show the very sur- 
prising ability of the nonadaptive protocol to demon- 
strate near Heisenberg-like scaling (SB ~ 1/T), even at 
low detection contrasts (/„ = 0.55). For low values of 
(Mk,F), the method exhibits a scaling which is worse 
than or comparable to the shot-noise limit. However, as 
these parameters are increased (Mk,F = 8 H> 12), the 
method exhibits a significant improvement and moves to- 
ward ~ 1/T scaling. In addition, this behavior occurs for 
a very large range of detection contrasts (f a > 0.55), and 
thus could be of great use to many other physical models 
besides the case of magnctomctry here, where one does 
not have perfect single-shot readout. 



V. CONCLUSION 

We have shown that unambiguous magnetic field 
strength estimation with sub-shot-noise scaling can be 
achieved by performing a time-varying sequence of mea- 
surements of the nuclear spin in a single N-V, using ei- 
ther adaptive or nonadaptive phase shifting of the system 
phase throughout the sequence. By numerically simulat- 
ing the protocol for both ideal and imperfect single-shot 
measurement conditions, we found the surprising result 
that the nonadaptive protocol, which is assumed to be 
far easier to realize experimentally, provides an estimate 
whose precision scales better than the adaptive proto- 
col. Another remarkable feature of our proposed method 
is that by sampling the phase accumulation at varying 
times the protocol yields an unambiguous estimate of 
the magnetic field. This means that one can determine 
an absolute estimate of the local magnetic field strength, 
even though the field may cause many Larmor rotations 
of the spin. The sub-shot-noise scaling of the estimate 
precision is only bounded below by the Heisenberg scal- 
ing (SB ~ 1/T), and the unavoidable dephasing of the 
spin. 

Recent experimental progress on coherent manipula- 
tions of the nuclear spin in the nitrogen- vacancy system 
in diamond allows one to perform SSM on the nuclear 
spin^ and this may soon lead to experimental demon- 
strations of our protocol with subsequent benefits for 
large scale nanomagnctic imaging. For the same total 
imaging time, the precision of our proposed estimation 
technique scales better than that of the standard esti- 
mation method. Moreover, our proposal is applicable to 
any other atomic-like physical systems, such as an elec- 
tron spin in silicon (SSM in that system has recently 
been reported^-), atomic vapors, or Bosc-Einstein con- 
densates. 
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